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Abstract 

Within the nonrelativistic QCD (NRQCD) factorization framework, we calculate the exclusive 
decay process rjb —* J/if^J/ifj to next-to-leading order in the strong coupling constant, while at 
leading order in the charm quark relative velocity. It is found that this new contribution to 
the amplitude is comparable in magnitude with the previously calculated relativistic correction 
piece, but differs by a phase about 90°. Including this new contribution will increase the previous 
prediction to B{r}\> — ► J/ipJ/ip) substantially, thus brightening the discovery potential of this clean 
hadronic decay channel of rjb in the forthcoming LHC experiment. 
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The quest for the rjb meson, the lowest-lying pseudoscalar bottomonium state, started 
shortly after the discovery of T in 1977. A thorough understanding of its properties, such 
as its mass and decay pattern, will be great benefit to strengthening our knowledge towards 
heavy quarkonium physics 

Tremendous efforts have been invested in various high energy collision experiments during 
the past thirty years to seek this elusive particle, but convincing evidence about its exis- 
tence has not emerged until very recently. After many futile searches, there finally comes 
out an exciting news from BABAR collaboration last month which reported the first 
unambiguous sighting of rjb with a 10 a significance, through the hindered magnetic dipole 
transition process T(3S) — > ^7. It is interesting to note the rather precisely measured 
hyperfine T(lS)-7] b mass splitting, 71.4±|;f(stat] ± 2.7(syst) MeV, can play a decisive role 
in discriminating different quarkonium models 

Aside from its mass, almost nothing is known about the dynamical properties of the rjb, 
e.g., its decay pattern. Owning to the weaker QCD coupling at bottom mass scale, together 
with the copious phase space opened up for innumerable decay channels, one in general 
expects that each rjb exclusive decay mode may only be allotted a rather small branching 
ratio. Indeed, a rough estimate of some 2-body and 3-body hadronic decay channels of rjb 
supports this expectation, and probably rjb disintegrates primarily into final states of high 
multiplicities yf. 

n 

Some "golden" modes have been proposed for hunting rjb, such as rjb — > J/ipJ /ip [4j and 
rjb — > J/ip'y @, E3]. Despite very clean signature thanks to the presence of J/if>, these decay 
modes unfortunately have rather suppressed branching ratios. Obviously, the e + e~ collision 
experiments like B factory are not well suited to pin down these decay modes due to the 
relatively low rjb yield. By contrast, high energy hadron collision experiments in general 
possess a much larger rjb production rate, which in turn allows for triggering these rare 
but very clean decay modes. For the decay rjb — > J/ipJ/ip, the original hope is that it 
should be discoverable in Tevatron Run 2 W] . Subsequent investigation reveals the situation 

n q 

may not be so optimistic [3|, [7] ■ In particular, an explicit NRQCD calculation predicts 



1 Note there recently appeared an estimate for this process based on a final-state interaction model, in 
which rjb first annihilates into DD , and the charm meson pair subsequently scatter into double J /if) Q . 
The corresponding prediction critically hinges on the branching ratio of 775 — > DD , which is poorly known 
at present. However, there is reasonable argument that suggests this ratio may be rather suppressed [^J. 
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B[?lb — > J/tpJ/ip] ~ 10~ 8 , which implies that notwithstanding the pessimistic prospect of 
observing this decay channel in Tevatron Run 2, the chance may still be bright in the 
forthcoming LHC program 3J]. 

The exclusive decay rjb — > J/ipJ/ip, apart from its phenomenological interest, is worthy 
of theoretical exploration in its own right. Since this hard exclusive process comprises 
entirely heavy quarkonium, it can serve as an ideal laboratory to test the applicability 



of NRQCD and perturbative QCD 



9j. Note this process is quite similar to the double 



charmonium production processes at B factory, e.g. e + e — > J/iprj c , J/ipJ/ip, which have 



stirred quite a few attention since the first Belle measurements were released in 2002 10 1 

nn r 



181, 



for an incomplete list about NRQCD-based investigation, see 

ro), [2^]). No doubt a systematic study of this hadronic decay process will offer complementary 
insight towards understanding the decay and production mechanisms of quarkonium. 

A peculiar feature of the decay rjb — > J/ipJ/ip is that 2 , at the lowest order in a s and 
charm quark relative velocity in J/if), v c , the amplitude vanishes. Therefore, to obtain a 
non- vanishing result, one must proceed to the higher orders either in a s or v c . In Ref. jsj, 
the leading relativistic correction to the tree-level process was explored. In this work, we 
make progress along the alternative route, i.e., we compute the 0(a s ) correction but stay at 
the zeroth order in v c . Since a s ~ the latter contribution is expected to be as important 
as the former one, so it will be definitely desirable to work out this missing piece. 

Before presenting our main result, it is instructive to recall some dynamical properties 
about this process. Let A and A stand for the helicities of two J / if) mesons in the rjb rest frame. 
According to the helicity selection rule and angular momentum conservation, the d ecay 
configuration (A, A) = (0, 0) exhibits the slowest asymptotic decrease with m b , B ~ l/wij [21 ]. 



However, it is impossible for rjb to decay into two longitudinally polarized J/ip because of 
conflict between parity and Lorentz invariance. The hadron helicity conservation, |A+A| = 0, 
can be violated either by the nonzero charm mass m c or by the transverse momentum of c 
inside J/ip, q±. For every unit of deviation of this rule, there is a further suppression factor 
of ml/ml or q\/m\. For the only physically allowed helicity configurations (A, A) = (±l, 
±1), the selection rule is violated by two units, so one expects B ~ 1/m 8 no matter the 



2 This decay can also proceed through rjb — > 7*7* — > J/tpJ/ip. Since this QED contribution is much less 
important than its QCD counterpart, it will be neglected in this work. 
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FIG. 1: Lowest order diagram, together with some representative next-to-leading order diagrams 
[from o) to e)] that contribute to % — ► J/ip J/ip. 



cause of violation is due to nonzero q± or nonzero m c . 

A noteworthy finding is that, at the lowest order in a s , the violation of the rule should be 
ascribable to the nonzero transverse momentum of c in J/ip, instead of the nonzero m c j^]. 
The corresponding branching ratio exhibits the following asymptotic scaling 

B v ,[ Vb ^ m±l) + m±l)} ~ alvt ( • ( 4) 2 ~ a]vf f^V , (1) 

\m b J \m z b J \m b J 

where the subscript emphasizes the corresponding contribution comes from the relativistic 
correction. The factor v\ affiliated with the would-be "leading-power" scaling, stems from 
the squares of the wave functions at the origin of two J/ip, since ipj/^p(0) oc (m c t> c ) 3//2 . In 
the last equation q± ~ m c v c has been assumed. 

One may naturally wonder how the situation becomes when one goes to the next-to- 
leading order (NLO) in a s . Some typical NLO Feynman diagrams contributing to r] b — > 
J/ipJ/ip are shown in Fig. [TJ Should it also be obligatory to retain finite relative momentum 
of c in order to obtain a non-vanishing result? To answer this question, it is helpful to 
draw some clue from an analogous hadronic decay process, T — > J/ipr] c , which also violates 
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the helicity selection rule. In this T decay process, the lowest-order diagrams start already 
at one-loop, i.e., at 0(a®), which in fact share the same topology as those represented by 
Fig. [Ib). There it is found that keeping a nonzero m c but neg lecting the relative momentum 
of c is sufficient to give rise to a non- vanishing result [22]. This may be reckoned as a 
persuasive hint that a similar pattern may occur in our case. A more direct evidence comes 
from inspecting the cut structure of Fig. [TJ It is easy to observe that all the diagrams in Fig. [TJ 
[except Fig. [T^)] permit a specific cut that has a clear physical meaning, i.e., a cut through 
two intermediate gluon lines that divides the full process into % — > gg and gg — > J/i/jJ/i/j. 
Since both sub-processes are finite at the zeroth order in v c , so does the corresponding 
portion of the imaginary part of the rjf, — > J/i/jJ/i/j amplitude which is obtained by stitching 
both on-shell sub-amplitudes together in accordance with the recipe of Cutkosky rule. Since 
a fraction of the imaginary part of the total amplitude is nonzero, there is no reason to 
expect the full imaginary part should vanish, thus we are reassured that at this order, the 
nonzero m c can adequately play the role of violating the helicity selection rule, and it suffices 
to set the relative velocity v c to zero in the calculation. As a result, we expect the NLO 
perturbative contribution scales as follows 

Ba.\m - w(±i) + w±i)] - <m (— y ■ <* (4V - « (— y > (2) 

\m b J \m£J \m b J 

where the subscript stresses the corresponding contribution comes from the 0(a s ) pertur- 
bative correction. As will be confirmed by the explicit calculation, this power-law scaling 
indeed holds besides logarithmical scaling violation effect. 

Having said this much, let us now proceed to the actual calculation. We first set up 
some notations. Let P,e (P,e) signify the momentum and polarization vector of each J/ip. 
Owning to parity and Lorentz invariance, the decay amplitude M. is constrained of the 
following structure: 

M = A P, P v 4(A) e£(A) . (3) 

Apparently, the only allowed helicity configurations indeed are (A, A) = (±1,±1). All the 
nontrivial dynamics is encoded in the reduced amplitude A, and our task is then to unravel 
its explicit form. 

In total 80 diagrams can be drawn for this process at 0{a & s ). A full treatment of these 
diagrams, at first sight, may seem to be a daunting task. Luckily, the calculation turns out 
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to be much simpler than one would imagine, thanks to some gratifying traits of this process. 
First, renormalization is not needed in this work. This originates from the fact that the 
amplitude at leading order in a s and v c vanishes. Since the contribution of the counterterm 
diagrams are proportional to the LO diagram, they simply vanish too. As a consequence, 
all the diagrams which contain primitive ultraviolet (UV) divergent subgraph must instead 
be free of the UV divergences, otherwise these unremoved divergences would signal the 
failure of the renormalization theory. A remarkable property of this process is that, these 
superficially UV-divergent diagrams are not only finite, in fact all of them vanish. Similarly, 
among the remaining manifestly UV-finite diagrams, a simplification further occurs, namely 
many diagrams with a symmetric topology also vanish, e.g., diagrams involving one 4-gluon 
vertex, diagrams involving two 3-gluon vertices, and diagrams containing a gluon ladder 
between the b and b lines. This can be mainly attributed to the specific Lorentz structure 
as indicated in ([3]), that is, in these diagrams, the Levi-Civita tensor that arises from the bb 
annihilation amplitude, is either contracted with some symmetric Lorentz indices, or there 
are no sufficient number of independent 4- vectors to contract with it. 

There finally survive 40 diagrams with non-vanishing contribution, some of which are 
shown in Fig. [T^,)-e). In fact, because of the interchanging symmetry between two J/ip, dia- 
grams of a given topology usually give rise to an identical answer, and practically only much 
fewer diagrams need to be calculated. We employ the automated Feynman Diagram Cal- 
culation package (FDC) to fulfill the analytic evaluation of the required one- loop diagrams. 
FDC is a powerful program based on the symbolic language Reduce, which is designed to 
automate the perturbative quantum field theory calculation in computer. FDC was initially 
developed by one of us long ago 23], and the function of automatic one-loop calculation 



has recently been realized by two of us 



24j. Recently, FDC has been successfully applied 



to several important quarkonium production processes to compute the NLO perturbative 



corrections 



19 



2Q|. 



In passing, we would like to comment on some specific issues in the calculation. One 
technical complication is the occurrences of six-point one- loop integrals in Fig. [Th) and b). 
FDC has implemented some systematic recursive algorithm to reduce higher-point integrals 
down to lower-point ones. In this work, it turns out that all the encountered 6-point integrals 
can be reduced to the 4-point integrals and lower ones. In addition, Passarino-Veltman 
algorithm 2Ja| has also been built into FDC to expedite the reduction of tensor integrals to 
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the scalar ones. 

It turns out that each individual diagram is infrared finite, at least in the Feynman gauge 
as we worked with. At a glance, one may feel that some care needs to be paid in han- 
dling Fig. [Tb), since the Coulomb divergences may potentially arise there. Hearteningly, a 
close examination shows that this diagram is in fact free of any singular l/v c poles. The 
absence of Coulomb singularity can be most transparently understood in the NRQCD fac- 
torization ansatz |9j . What we can calculate within perturbation theory is the so-called hard 
(short- distance) coefficient. All the long-distance singularities encountered in computing the 
one-loop QCD diagrams, must be exactly reproduced in the respective one-loop NRQCD 
diagrams with a potential gluon attached between the outgoing c-c lines that forms a J/ip, 
so they ultimately cancel out in the intended hard coefficient via the matching procedure. 
However, the corresponding diagrams in the NRQCD sector trivially vanish in our case, 
again because the % — > J/ipJ/ip process vanishes at LO in a s and v c , therefore no any net 
long-distance divergence, such as Coulomb divergence, should be expected to appear in the 
diagrams in the QCD sector. 

Finally it may be worthwhile to compare our process with T — > J / ipr] c . As stated earlier, 
the latter process is described by a small subset of Fig. [TJ the class of diagrams exemplified by 
Fig. [lfe), in which three "Abelian" gluons connect between b quark line and the charm quark 
lines. Since T has negative C parity, charge conjugation invariance demands that three 
gluons must form a totally symmetric color-singlet state, therefore only the fully symmetric 



part of the color factor d abc d abc can contribute in this T decay process 22j. On the contrary, 
due to the positive C parity of rjb, in our case the three gluons must instead be totally 
antisymmetric in color space, so the f abc f abc piece finally survive in the color factor. We 
have explicitly checked this is indeed the case. 

After analytic manipulation of all the diagrams by FDC, together with some simplifying 
by hand, we can express the reduced amplitude as following: 

512^ 



9^^(0)^(0)/ -f , (4) 



where the relativistic normalization is used for the quarkonium states, and |P| = yj ml — 4m 2 
is the magnitude of J/tp momentum. The complex- valued / function, which has encom- 
passed the loop contribution, is deliberately chosen as dimensionless, therefore it can depend 
on m& and m c only through their dimensionless ratio m 2 c jm\. Furthermore, the / function 
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FIG. 2: The solid lines represent the real and imaginary parts of /(£), the dashed lines represent 
their respective asymptotic expressions as given in ([6]) and 0. The vertical line marks the place 
£ = m^/ml = 0.102, which is adopted in the phenomenological analysis. 



is normalized such that it varies with m 2 c / only logarithmically, which is slower than any 
power law scaling. After some quick algebra, one can verify that Eq. (Tj0) is in conformity to 
B ag ~ 1/mf up to logarithmical modification. 

We display the shape of the / function in Figure [2J The analytic expression for the real 
part of the / function is too lengthy to be reproduced here. However, the analytical result 
for the imaginary part is quite compact, 



Im/(0 



7T 



1 



P _ 63 + 86£ 2 
2£~ 2(2 + + P 



23 

T 



- 14M tanh" 1 (3 



8 U(2 + 

+ (53 - 140 tanh- 1 (0 - (1 - 20 ln(2 - 40 - (1 + 20 ln(4£) } 



(5) 



where (3 = ^1 - 4£. 

It is instructive to deduce the asymptotic behavior of / in the small £ limit. After some 
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efforts, we obtain 



lm/ asym (0 



R'6/ aS y m (^) 




(6) 



(7) 



As can be clearly visualized in Fig. [2, as £ — >• 0, the real and imaginary parts of the / function 
exhibit ln 2 £ and ln£ dependence, respectively. This seems to be a generic trait of exclusive 
bottomonium decays to double charmonium, which is also seen in the analogous / function 
in the T — > J/V^c process [22]. It is interesting to note that, up to the physically relevant 
point £ = ml/ml = 0.102, which is obtained by substituting m c = 1.5 GeV and nib = 4.7 
GeV, the asymptotic expressions in (jBJ) and ((Zj) still constitute decent approximations to the 
true expressions, where the agreement for the real part is especially satisfactory, and the 
agreement for the imaginary part is at 70% level. This may suggest that, if one is willing to 
carry out this calculation by hand yet without loss of accuracy, the most efficient way is to 
first expand the integrand in vr? c jm^ prior to performing the loop integrals. 

It is interesting to note that, provided that £ is not overly small, say, £ > 3.8 x 10 -4 , 
then — Im/ is always bigger than |Re/|. In particular, at the phenomeno logically relevant 
point £ = 0.102, we find the corresponding / = 0.216 — 2.30i = 2.31 e _l84 ' 6 °, and the real 
part seems to be practically negligible. 

Since the NLO perturbative amplitude is dominated by its imaginary part, one may 
like to speculate which type of cuts might constitute the most important contribution. It 
might sound attractive to inquire whether that specific cut as mentioned earlier, the one 
acting upon two intermediate gluon lines that partitions the full process into T]b — > gg and 
gg — ► J/ijjJ/ijj, would dominate over all other possible cuts. If this were the case, we 
would be happy because this process then could be endowed with a "plausible" physical 
picture - that it proceeds first through r]b annihilation into two gluons, and these two gluons 
subsequently rescatter into double J/ip. Let us examine this conjecture concretely. With 
the aid of the Cutkosky rule, it is straightforward to deduce this particular fraction of the 
imaginary part of the / function, which we call Im/ 2 g : 
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Note the occurrence of the l/£ term, which can be traced to the cut contribution from 
Fig. [T^l) , renders this function scale with quark masses as ml /ml, so that the true scaling 
B as ~ I /mf is violated. This is unacceptable and clearly indicates approximating the full 
imaginary part by this specific fraction is unjustified. Numerically for £ = 0.1, we obtain 
Im/29 — 115, which is far off the true value Im/ = —2.30. A careful investigation reveals 
that each diagram in Fig. [1] [except Fig. [Ife)] possesses a rich cut structure, for example, 
Fig.[[k,), b), d) allow for 6, 5, and 4 different kind of cuts, respectively. Among all the possible 
cuts, those which pass through one gluon line and another nonadjacent quark line generally 
lead to more involved expressions in (jSJ). For a given diagram, it seems not profitable to 
single out one specific cut to mimic its full imaginary part. For example, only after the 
remaining cuts are included, the imaginary part of Fig. dh) then has a sensible answer, free 
of the dangerous l/£ singularity. 

We are now very close to the intended formula. For completeness, we still need the 
expression of the reduced amplitude that arises from the tree-level relativistic correction jj] 

, 512v / 6tt 2 oil m c , . 2 . , n . 

M = 1 13/2 ^(0) Vj/^(0) (v 2 ) m ■ (9) 

27 m b 

Roughly speaking, (v 2 )j/ip is a quantity related to the second derivative of the wave function 
at the origin for J/ip, which governs the size of relativistic corrections. In fact, this factor 
admits a rigorous definition as a ratio of NRQCD matrix elements 

where e denotes the polarization vector of J/ip at its rest frame, ip and \ represent Pauli 
spinor fields in NRQCD, D is the spatial covariant derivative. The matrix elements appear- 
ing in the above ratio are understood to be the renormalized ones. It is worth stressing that 
the sign of {v 2 )j/^ needs not necessarily to be positive. 
Plugging equations PJ and © in the formula 



IP 



3 



in 



we then obtain the desired partial width. Note the cubic power of momentum reflects that 
two J/ip are in relative P-wave state. This formula has already taken into account the sum 
over transverse polarizations of two J/ip as well as their indistinguishability. It may be more 
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convenient to present an explicit expression for the branching ratio, where the wave function 
at the origin of r)b drops out, 

•> 13 7r 2 a? ml , A _ A 4m^ 3/2 



B[ Vb -> J/> J/V>] = iC — ^ -f ^m(O) 



a, 4m? \ 1 „ / m 2 
- \ mi ) \m 2 



(12) 



In deriving this, we have approximated the total width of rjb by its hadronic width: 



r to t h>] « r[?7 6 -> light hadrons] = — -f < (0) , (13) 



3m 2 



where the factor K gg encodes the magnitude of NLO perturbative correction to rjb — > gg 26] 



'53 3l7r 2 8nA a s (2m h ) 
K 99 = 1 + I — 



2 24 9 J 7i 
rif denotes the number of active light flavors. One usually takes rif = 4, by treating charm 
quark also as a massless parton. For a s {2mb) = 0.18, we get K gg = 1.58. 

Equation (|T2l constitutes the key formula of this work. From this equation one can 
readily confirm the asymptotic behaviors anticipated in (CQ) and (j2j). 

Let us now explore the numerical outcome of (TI21) . The input parameters are mb, m c , 
a s , ipj/^(0) and (v 2 )j/^, respectively. We take m& = M Vb /2 ~ 4.7 GeV and m c = 1.5 GeV. 
Obviously, our predictions are quite sensitive to the value assigned for the strong coupling 
constant. To account for the affiliated ambiguity, we slide the renormalization scale [i from 
2rrib to 2m c , corresponding to varying a s between 0.18 and 0.26. The wave function at the 
origin of J/ip can be extracted from its electric width: 

r m e+e-] = ^f^(O) (l - , (14) 

where the first-order perturbative correction has been included. Using the measured electric 
width 5.55 keV [27], we obtain ipj/^O) = 0.263 GeV 3/2 . 

Among various input parameters, the least known about is (v 2 )j/^. Nevertheless, there is 
a useful relation, first derived by Gremm and Kapustin by employing the equation of motion 
of NRQCD [28], which expresses this quantity in terms of the charm quark pole mass and 
the J /if) mass. To our purpose this relation reads [11] 



M 2 Jf -4m 2 cpole 

\ v )m ~ A m 2 • ( 15 ) 
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FIG. 3: B[r]b — ► J/ifiJ/i/j] as a function of (v 2 )j/^ (left panel) and a s (right panel). The solid 
lines represent the full prediction from (|12p and the dashed lines are obtained by keeping only the 
relativistic correction piece. In the left panel, three different pairs of curves (each pair tinted by the 
same color) are obtained by fixing the renormalization scale [i to be 2m&, m;, and 2m c respectively, 
so that the corresponding strong coupling constant a s is 0.18, 0.22 and 0.26. 



Owing to the intrinsic ambiguity associated with the charm quark pole mass, a wide spec- 
trum of estimates for m cpo i e is scattered in literature. Ref. pjj takes m cpole = 1.4 GeV in 
the analysis of double charmonium production at B factory, which leads to (v 2 )j/^ = 0.22 
from (tl5l). Note this value is compatible with (v 2 )jm = 0.225+o;o88, 

which is obtained from 

a recent Cornell potential-model based analysis fl Q- A positive (v 2 ) j/j/, is helpful to 
alleviate the discrepancy between the LO NRQCD prediction and the measured J/ip + i] c 
production rate at Belle. Notwithstanding the phenomenological inclination, however, one 
may be alert that a rather different value might not be excluded on theoretical ground. For 
example, if m cpo i e = 1.75 ± 0.15 GeV is assumed, which stems from a recent QCD moment 
sum rule analysis [3l| . one would instead obtain (v 2 )j/^ = —0.22 ± 0.15, whose sign has 
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even changed. Because of the impossibility to nail down m cpo i e to an accuracy better than 
Aqcd, it seems safe to take a conservative attitude, by assigning a rather large uncertainty 
to (v 2 )j/^. In this work, we will assume —0.3 < (v 2 )j/^ < 0.3, including the possibility that 
it may even vanish. 

In Figure[3]we illustrate how the branching ratio of rjf, — *> J/ip J/ip varies with (v 2 )j/^ and 
a s . In the left panel of Fig. [3j one can clearly see that inclusion of the QCD perturbative 
correction will significantly enhance the previous prediction that only implements relativistic 
correction, especially when \(v 2 )j/^\ is small. The approximately parabolic dependence of 
B on (v 2 )j/,p, reflects the fact that the NLO perturbative amplitude, A aa , is dominated by 
its imaginary part, therefore it is nearly orthogonal to the relativistic correction amplitude 
A v 2, so their contributions can be added almost incoherently. To see the interference pattern 
lucidly, we make a numerical presentation to the terms inside the absolute value sign in ffT2l , 

(v 2 )j/^ + 1.24e- i84 - 6 °a s = ((v 2 ) J/f + 0.12a s ) - i 1.24a,. 

It is interesting to note, besides the approximate 90° phase difference, the (imaginary part 
of) NLO perturbative contribution is comparable in magnitude with the relativistic correc- 
tion piece. The residual interference between the real part of the perturbative correction 
amplitude and the relativistic correction piece has modest effect, which can be either con- 
structive or destructive depending on the sign of {v 2 )j/^. In the right panel of Fig. [3l one 
can further observe that, despite the acute sensitivity to the choice of the renormalization 
scale, the complete result is considerably larger than each individual contribution in the 
considered range of a s . 

As can be inferred from Figure [3j our final estimate of the branching ratio is 

B[r} b -> J /if) J/ip] = (2.1 - 18.6) x 10~ 8 . (16) 

Not surprisingly, the uncertainty is rather large. It is enlightening to compare our result 
with the previous prediction that includes only the relativistic correction contribution j^]: 

B V 2 [r) b -> J/ip J/ip) = (0.5 - 6.6) x 10~ 8 . (17) 

Though the strategies of estimating the error differ in both work, it is nevertheless un- 
ambiguous to see the substantial effects of the NLO perturbative correction. It is worth 
mentioning that the lowest bound of our prediction is already as large as 2 x 10~ 8 , which 
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can be entirely attributed to the NLO perturbative contribution. In contrast, the prediction 
including only the relativistic correction simply vanishes as (v 2 )j/^ approaches zero. 

Finally we present an updated estimate about the discovery potential of this decay mode 
in the LHC experiments. In the hadron collision experiment J/ip can be cleanly recon- 
structed via its muonic decay mode. Multiplying (I16I) by the branching ratios of 6% for 
each of the decay J/ip — > we obtain B[r] b — > J/ipJ/ip — > 4/x] rs (0.7 — 6.7) x 10~ 10 . 

Assuming the rjb production cross section at LHC to be about 15 /xb |3j, one then finds the 
cross section for the 4 ji events to be about 1-10 fb. For a 300 fb _1 data, which is expected 
to be gleaned in one year run at LHC design luminosity, the number of produced events may 
reach 300-3000. The product of acceptance and efficiency for detecting J/i[) decay to ji + ji~ 
is estimated to be e « 0.1[4|, which is perhaps a conservative estimate for LHC. Multiplying 
the number of the produced events by e 2 , we expect between 3 and 30 observed events per 
year. From this rough estimate, we are tempted to conclude that, the prospect of observing 
rjb at LHC through the 4/i mode looks promising. However, we should be fully aware that it 
will be a highly challenging task for experimentalists to single out the relatively few signal 
events from the presumedly abundant continuum background. 

To summarize, in this work we have studied the exclusive decay process rjb — > J/ipJ/ip to 
NLO in a s while at LO in v c . We found that this new contribution to the amplitude is com- 
parable in magnitude with the previously calculated tree-level relativistic correction piece, 
but differs by a phase about 90°. Including this new piece of contribution will substantially 
enhance the previous estimate of B(r]b — > J/ipJ/ip), thus the observation prospect of this 
clean rjb hadronic decay mode in the forthcoming LHC experiment becomes brighter. 
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